Random spreading phenomena in annealed small world networks 
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We study the simple random walk dynamics on an annealed version of a Small- World Network 
(SWN) consisting of A'' nodes. This is done by calculating the mean number of distinct sites visited 
S(ri) and the return probability Pooit) as a function of the time t. S{t) is a key quantity both from 
the statistical physics point of view and especially for characterizing the efficiency of the network 
connectedness. Our results for this quantity shows features similar to the SWN with quenched 
disorder, but with a crossover time that goes inversely proportianal to the probability p of making 
a long range jump instead of being proportional to as in quenched case. We have also carried 
out simulations on a modified annealed model where the crossover time goes as due to specific 
time dependent transition probabilities and we present an approximate self-consistent solution to it. 
PACSnumbers: 05.40.-a, 05.50.+q, 87.18. Sn 



I. INTRODUCTION 



During the last few years an overwhelming amount of 
evidence has been accumulated about diverse networks 
showing "small world" properties. This means that on 
the average an arbitrarily selected node can be reached 
from another node in very few steps despite the fact that 
only relatively small number of links are present [Q . Such 
networks include social nets , the internet Q , the www 
of linked documents scientific cooperation relations 
etc. Watts and Strogatz were the first to suggested 
a simple mathematical model, which refiects the small 
world phenomenon: Take a regular lattice and introduce 
(few) random links between any pairs of sites Q. This 
Small World Network (SWN) model interpolates between 
regular graphs or lattices and the so called Erdos-Renyi 
random graphs [Q. 

Real networks are commonly characterized by a number 
of parameters but perhaps the most important measure is 
the average distance between their nodes. It has turned 
out that there is a rich family of models with the small 
world property but they differ in many other respects. 
For example, the degree distribution of the nodes is Pois- 
sonian for the SWN while many real networks are often 
scale free (with a power law decay in the degree distri- 
bution) for which preferential growth models have been 
introduced §. Thus SWN's are very interesting graphs 
not only because of these properties but also because they 
are simple with nontrivial behavior but provide possibili- 
ties for exact solutions |9j and they seem to have a direct 
application in polymer physics |10|| . 
I addition to static structural properties of networks there 
is ever growing interest in the dynamics processes in net- 
works. In spite of this, relatively few results have sofar 
been published on the dynamics of small world models. It 
is expected that the underlying network topology should 



have a major impact on practically any phenomenon tak- 
ing place on it. This is supported by the recent results on 
the the spectral density of the adjacency matrix of small 
world models, which show that these graphs produce a 
dramatic deviation from the semi-circle law of random 
graphs |ll|. 

It is our belief that spreading phenomena are perhaps 
the most direct examples of dynamical processes reflect- 
ing the small world properties. In direct spreading of 
e.g. a disease the nodes of the graph get infected by the 
rule that infection propagates each time step to all un- 
infected neighbors of already infected sites iQ. Then 
the simplest example of non-trivial dynamics could be 
that of a diffusing particle on the SWN. This in turn 
is related to the intensively studied process of random 
walks in random environments p3[ |. Recently some re- 
lated papers have been published, for example the study 
of spectral properties of the Laplacian on the SWN . 
While some numerical and analytic results were given in 
p6| for the spreading phenomenon being characterized by 
the average access time to the sites of the system, Jasch 
and Blumen |l4| published simulation results for spread- 
ing on SWN using random walk dynamics with the main 
quantity of interest being the average number of distinct 
sites vistited at a given time. 

In all these studies a simple version of the SWN was con- 
sidered: A one-dimensional ring consisting of N nodes 
with fc-neighbor interactions was assumed where, in ad- 
dition, Np new links were introduced between arbitrar- 
ily chosen, but not yet connected nodes. The result by 
Jasch and Blumen |14[| is intuitively very clear: whenever 
there is however small ratio of long range links present in 
the graph, S{t) changes its asymptotic behavior from the 
~ y/i relationship characteristic for the one-dimensional 
case io ^ t representing the long range interconnected 
graph. In a crossover formula was presented by re- 
lating the probability of long range links p to a character- 
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istic time t*, where the crossover sets in =const). 
Again, one would intuitively expect that a = 2 and that 
the numerically found discrepancies were occuring most 
probably due to limited system sizes. In fact we found 
psf in our numerical simulations that a ~ 2. This turned 
out to be true not only for SYn) but also for the so called 
return probability Pqo 1T5|JT7[ |. 

In their original form the SWN is defined with quenched 
disorder represented by the long range links. However, it 
is interesting to ask to what extent the small world prop- 
erties of the dynamics depend on the quenched character 
of disorder. In order to find some answers we are in 
this paper setting to study the annealed version of the 
dynamics on SWN's, where the links are not frozen in 
but the walker has at any site the choice to make either 
a long range jump (with probability p) or to continue 
the random walk in the regular part of the lattice (with 
probability (1 — p))- Such a model was studied in this 
context by Pandit and Amritkar where the average 
access time was calculated. Here we will concentrate to 
determine through analytical theory and simulations the 
average number of distinct sites visited and the return 
probability. 

The paper is organized as follows: In the next section we 
introduce the model and present the solution. In Sec- 
tion 3 we discuss a modified version of the model where 
the crossover time scales as p~^ as in the quenched case. 
Finally we present a discussion and a summary. 



W = (1 -p)W(^) +pW 



(3) 



Here (S) and (L) refer to short and long range jumps, 
respectively. The zeroth row of W^'^^ reads as follows 



k times N-2k-l times k times 
A similar equation holds for the W'^h 
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,0,...,0). 



fc+1 times N-2k-l times k times 



(5) 

The i*^ rows are then obtained by cyclically shifting the 
0*'* rows to the right. All W and T matrices have the 
Toeplitz form, i.e., Ty depends only on the difference 
{i — j). Therefore, the right hand side of Eq. (1) is a 
convolution which leads after spatial Fourier transform 
to the following form 



With the initial condition 

the formal solution is as follows 



(6) 



(7) 



II. THE BASIC MODEL AND ITS SOLUTION 



P,(i) = exp[ / {Wgiu) - l)du]. 



(8) 



In order to visualize the model let us consider a ring 
with first and second neighbor connections (fc ~ 2), the 
sites of which being numbered from to A^ — 1. In this 
system we let a random walker to start at site at time 
t = 0. Before the first step sets in the four links of 
site are rewired by moving the other endpoint to any 
of the other sites with probability p. Then the walker 
moves to one of the sites newly linked to its position 
with equal probability. Let us call this site i. After the 
walker arrives the site i the original lattice is restored, 
the four links for i are rewired as before, etc. In this way 
the small world property is maintained but we have got 
rid of the problems of quenched randomness. 
The movement of the random walk is governed by the 
simple master equation: 



^tP^{t) 



T^,P^{t) 



(1) 



where the continuum time limit has been applied. Here 
Pi (t) is the probability that the walker is at site i at time 
t and 



(2) 



with Wij being the transition matrix, which has the fol- 
lowing form: 



This solution can be easily evaluated for the matrix W 
given in Eqs (3)- (5). 



A. Return probability and spreading rate 

Let Fij {t) denote the probability of the random walker 
visiting site j at time t having started from site i. Then 
we can write 



Jo 



(t — u)du, 



(9) 



where Pij (t) is the probability for the random walker to 
move from i to j during the time t. From this we get 
through the Laplace transform the following equation 



(10) 



Now let us take s{t) as the probability of observing a new 
site, or as the spreading rate at time t when the random 
walker started from site 0: 



s{t) 



E 

i=0 



Fo^it) 



(11) 
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By taking the return probabilities Puit) to be the same 
for all i the equation reads as follows 

1 1 

As stated before the quantity of interest is the average 
number of distinct sites visited, S(t), which is expressed 
as 

S{t) = / s{u)du. (13) 
Jo 

Then by using Eq. (|2|) we obtain S{t) by the inverse 
Laplace transform of the function s{z)/z. It is noted 
that S{t) could physically be intepreted as the amount 
of spreading the random walker has covered at time t. 



B. Results 



In order to obtain the return probability and the mean 
number of visited sites as a function of the probability 
p we have solved numerically the set of equations (||) 
and (|l^) given above. For the purpose of checking the 
validity of these solutions we compare them with discrete 
time simulations both for Pao{t) and S{t). As depicted 
in Figs |l| and || the agreement for both these quantities 
is excellent, apart from the short time behavior, which is 
a consequence of the difference between the discrete and 
continuum time versions of the system. 
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FIG. 1. The simulation results for the return probabilities 
forp= 10"*, 10"^■^ 10"^ 10"^ '^ (uneven line). The smooth 
curve shows the result of the analytical theory. 
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FIG. 2. The scaled return probabilities against scaled time 
for p = 10"'', 10"^ ■^ 10-^ 10-2 ■^ 
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FIG. 3. The simulation results for the mean number of dis- 
tinct sites visited S{t) (i.e. the amount of spreading) for 
p = lO""*, IQ-^-^ IQ-^, 10-2-*. These curves start from 
5(0) — 1. Analytical results, which start at 5(0) — 0, are 
also shown. 

Fig. |l| shows the return probability Pqq as a function 
of time. It is clearly evident from this figure that the 
simulation results agree very well with analytical calcu- 
lations. In Fig. H we present a scaling plot of this data, 
i.e. (Poo — l/A^)'\/i vs. scaled time pt, where 1/A^ repre- 
sents the long time limit. This scaling obeys the following 
functional form: 

Poo(t)=t~^/V(ip") (14) 

where (/3 is a universal scaling function with the the prop- 
erty that (p{x) oc const for x ^ 1 and it decays rapidly 
for a; > 1. The exponent a is unity (see Fig. H) which 
is easily understood, because the short time characteris- 
tic for the one-dimensional case is valid until the walker 
makes a long jump. This takes the time 1/p, on the 
average, which goes as a characteristic time into the ar- 
gument of the function /. This is in contrast with the 
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behavior of the quenched model, where in order to make 
a long jump, the distance proportional to \/p has to 
be abbridged by diffusive motion on a one-dimensional 
topology of the system leading to a characteristic time 
behavior of ~ {l/pf @,|l5|. 




pt 

FIG. 4. The scaled spreading {S{t)/y/i) against the scaled 
time (pt) forp= 10"*, 10-^■^ 10-^ IQ-^ ^ 

In Fig. ^ we show the simulation and analytical results 
for the mean number of distinct sites visited, i.e. S(t) as 
a function of time. Again, the good agreement between 
the simulation results and theoretical analysis is clear. 
The scaling form corresponding to (|l^) would now read 
as 

S{t) = t^'^fitp^) (15) 

where the universal scaling function / is of the following 

form 

J., . [const for x <C 1 

^ for x»l (1^) 

Here the short time behavior is 5 t^/'^ while for long 
times it becomes proportional to t. As seen from the 
scaling plot of Fig. |4| the same a = \ applies to this case 
as for the return probability. 



III. THE SELF-CONSISTENT MODEL 

Next we ask the question whether it is possible to modify 
the annealed model in a way that the crossover is shifted 
such that a = 2 would be obtained. This behavior is 
achieved by assuming a time dependent transition prob- 
ability, as discussed below. 

Since the scaling of the transition occurs in the quenched 
system later (as ~ p^t) we replaced the multiplier p of 
W^^^ in Eq. (||) with ps{t) to simulate the situation 
where the walker has a probability of making a long range 
leap only when visiting a previously unseen site: 



W = (1 -p)W(^^ -(-ps(i)W(^). (17) 

The corresponding master equation cannot be solved ex- 
plicitly but we can still estimate the solution with arbi- 
trary accuracy by numerical iteration scheme. In Fig. 
^ we show the resulting time dependent behavior of the 
random walk spreading for the self-consistent model and 
simulated quenched system. Apart from the short times 
the agreement between these results seem to be quite 
good. 
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FIG. 5. The results for the spreading vs. time of the 
self-consistent model obtained from the analytical theory 
(solid line) and from the quenched simulations (dashed line) 
for p = 10"'', 10"^ ■^ 10-^ 10-2 ^ 

In Fig. H we show the scaling plot with a = 2 for the 
self-consistent model. Apart from early times the scaling 
seems to hold. 
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FIG. 6. The scaled spreading of the self-consistent model 
S{t)/Vt against scaled time p^t for p = 10"*, lO"^'^, 10"^, 
10"^-^ 

Hence the numerical results justify the choice of Eq. ( p^ ) 
which reflects the fact that in the quenched model one 
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dimensional random walk has to be carried out between [17] A. Scala, L. A.N. Amaral and M. Barthelemy, cond- 
two long jumps. mat/0004380 

[18] R. Monasson, Eur. Phys. J. B 12, 555-567 (1999) 



IV. SUMMARY 



We have shown that the annealed random walk model 
with rare long jumps reflects some aspects of the 
quenched SWN's. In the simplest case, with time inde- 
pendent transition probabilities, the model can be solved. 
However, as expected, only qualitative agreement be- 
tween the quenched and the annealed models can be ob- 
served. With properly chosen time dependent transition 
probabilities we obtain even the proper exponent a — 1. 
We think that it is quite interesting that the behavior of 
the spreading can be estimated by an annealed model. 
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